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ABSTRACT 

Extended halo tidal streams from disrupting Milky Way satellites offer new opportunities for gauging 
fundamental Galactic parameters without challenging observations of the Galactic center. In the 
roughly spherical Galactic potential tidal debris from a satellite system is largely confined to a single 
plane containing the Galactic center, so accurate distances to stars in the tidal stream can be used to 
gauge the Galactic center distance, Rq, given reasonable projection of the stream orbital pole on the 
Xqc axis. Alternatively, a tidal stream with orbital pole near the Ygc axis, like the Sagittarius stream, 
can be used to derive the speed of the Local Standard of Rest (Glsr)- Modest improvements in current 
astromctric catalogues might allow this measurement to be made, but NASA's Space Interferomctry 
Mission (SIM PlanctQuest) can definitively obtain both Rq and Glsr using tidal streams. 
Subject headings: Milky Way: structure - Milky Way: dynamics - Sagittarius dwarf galaxy 



1. DISTANCE TO THE GALACTIC CENTER 

With the assumption that globular clusters trace the 
general shape and extent of the Milky Way (MW), Shap- 
ley (1918) first showed how they can be used to estimate 
the distance (Ro) to the Galactic center (GC), expected 
to lie at the center of the cluster distribution. Though 
Shapley's first execution of this experiment exaggerated 
Ro due to cluster distance scale problems, the overall 
scheme of mapping an extended distribution of Galac- 
tic tracer objects to determine the location of its center 
remains a valid, if traditionally underutilized, strategy. 

The globular cluster sample is relatively small and con- 
centrated to the GC, where dust effects introduce large 
distance uncertainties and a likely still incomplete and 
lop-sided cluster census. Population II tracers like RR 
Lyrae, blue horizontal branch (BHB) or giant stars are 
much more plentiful outside of the MW bulge. Unfortu- 
nately, the current census for these tracers is even more 
incomplete than for globulars. Though this situation 
may be remedied by currently planned wide angle sur- 
veys, several inherent problems remain with exploitation 
of these tracers as GC benchmarks. As with the clusters, 
the MW Zone of Avoidance (ZA) will always result in 
biased sample distributions and potential i?o underesti- 
mates — exacerbated if surveys do not reach the far side 
of the MW. Even more challenging is that the global dis- 
tributions of halo stars are far from dynamically mixed: 
Recent surveys of the above tracers reveal a halo streaked 
with substructure, likely originating as satellite disrup- 
tion debris (e.g., Vivas et al. 2001, Newberg et al. 2002, 
Majewski 2004) and eroding simple halo axisymmetry. 

The very existence of numerous tidal streams moti- 
vates the present contribution. Individual tidal streams 
actually possess a relatively simple spatial configuration. 
Within spherical potentials, tidal debris arms from a dis- 
rupting satellite will lie along the satellite orbital plane, 
which contains the GC. A sufficiently extended tidal de- 
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bris arc defines that plane, which intersects the MW Xqc 
axis at the GC. This simpler, almost two-dimensional 
geometry of tidal stream arcs removes the need for sam- 
ple completeness: In principle, Rq should be derivable 
from the (I, b, distance) distribution of only a large enough 
sample of tidal stream stars to define their orbital plane. 3 
In reality, non-spherical potentials precess tidal 
streams. Fortunately, this is a relatively small effect in 
the MW, at least for Rgc's of tens of kiloparsecs. John- 
ston et al. (2005; "J05") showed the Sagittarius (Sgr) 
tidal stream precession is sufficiently small to conclude 
the MW potential is only slightly oblate within the Sgr 
orbit (peri:apo-Galacticon of 13:57 kpc). Moreover, as 
pointed out also by Helmi (2004), that part of the Sgr 
trailing arm arcing across the southern MW hemisphere 
(see Majewski et al. 2003, "MSWO" ) is so dynamically 
young that it hasn't had time to precess (see Fig. 5 of 
J05). Unfortunately, as noted by MSWO, Sgr is in almost 
the worst possible orientation to undertake the proposed 
i?o-gauging: With a virtually negligible angle between 
the Sgr orbital plane and Xqc axis, small errors in the 
definition of the orbital plane (due to small residual pre- 
cession, and the finite width of the debris plane) lead to 
substantial uncertainties in derivation of Rq. 

The ideal tidal debris configuration for estimating Rq 
has a pole closer to the Xqc axis. Given the pace of 
discovery, such a stream may soon be found. Based on 
the nearly polar orientation of the HI Magellanic Stream 
and the typically measured proper motions (yu's) for the 
Magellanic Clouds (Gardiner & Noguchi 1996, van der 
Marel et al. 2002 and references therein), it is clear that a 
stellar counterpart to the Magellanic Stream would have 
almost the perfect orientation for gauging Rq. 

Systematic errors in a tracer distance scale translate to 
estimates of Rq. However, because streams contain dif- 
ferent stellar types (e.g., giant stars, RR Lyrae, BHB), 
uncertainties from photometric/spectroscopic parallaxes 
can be cross-checked. In most cases, reddening and 
crowding effects can be of negligible concern. Alterna- 
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tively, with NASA's Space Interferometry Mission (SIM), 
direct trigonometric parallaxes will be well within reach: 
For a putative Magellanic stellar stream orbiting at ~ 50 
kpc radius, the ~ 10-20 /xas parallaxes are well above 
the SIM wide-angle astrometric accuracy goal of 4/ias, 
assuming K giant star tracers (V ~ 18). 

As a test of what might be achieved, we ran N-body 
simulations of different mass satellites disrupting for 5 
or 10 Gyr (whatever was needed to produce > 270°-long 
tails) in the Galactic potential that best fitted the Sgr 
debris stream in Law et al. (2005; "L05" hereafter). The 
orbit was constrained to match the current position, ra- 
dial velocity (RV) and \i (Gardiner & Noguchi 1996) of 
the Small Magellanic Cloud (SMC), with orbital pole 
(l,b) — (196, —5)°. 4 All other model parameters were 
similar to those in L05. Each simulation was "observed" 
outside a \b\ > 15° ZA, with - 10 3 , 10 4 and 10 5 tracer 
stars (apportioned with ~ 90% of these in the satellite 
core and ~ 10% in the tidal tails), and with 0, 10 and 
20% random Gaussian distance errors imposed. The sim- 
plest (though not best!) analysis of these data is simply 
to fit a plane and measure its intersection with the Xqc 
axis. With this Oth-order method, even for large sam- 
ples of stars in dynamically cold streams (i.e., not that 
from a 10 10 M Q progenitor) and no distance errors, rel- 
atively large (< 7%) systematic errors in Rq can remain 
(Fig. 1) because plane-fitting does not account for the 
residual precessional twisting of the debris arms. The 
direction of precession (determined by the direction of 
the stream angular momentum vector) drives the sense 
of the imposed systemic Rq error (i.e. closer or farther), 
and random distance errors add additional uncertainties 
depending on details of the stream orientation relative to 
the ZA. A better, now proven method (e.g. L05) is to use 
N-body modeling to reconstruct a given stellar stream; 
such modeling can precisely account not only for pre- 
cession but also for stream dispersion and other higher 
order uncertainties, which would permit a more accurate 
identification of the center of the MW potential for an 
appropriately oriented stream. 

Recent measurements of stellar motions around Sgr A* 
have led to dynamical parallaxes good to 5% (7.94 ±0.42 
kpc; Eisenhauer et al. 2003), a measurement sure to im- 
prove with longer Sgr A* field monitoring campaigns. 
Few percent quality trigonometric parallaxes of stars 
near the GC will be measured as part of a SIM Key 
Project. In either method, the target stars are reason- 
ably expected to lie at the assumed center of the MW 
potential. The proposed use of tidal streams to measure 
i?o will provide an interesting test of this hypothesis, 
since tidal streams orbit the true dynamical center of the 
integrated potential over tens of kiloparsec scales. A com- 
parison of this center to the Sgr A* distance could reveal 
whether the MW may be a lop-sided spiral (e.g., Bald- 
win et al. 1980, Richter & Sancisi 1994, Rix & Zaritsky 
1995). Such lop-sidedness can, in fact, be induced by 
mergers of large satellites (Walker, Mihos & Hernquist 
1997). In principle, three well- measured tidal streams 
can verify whether the GC lies along (I, b) — (0, 0), since 
the true GC should lie at a mutual intersection of the 

4 We do not model the possibly complex interaction between 
the Small and Large Magellanic Clouds since we arc interested in 
testing a hypothetical stream with desirable properties. 



three corresponding stream orbital planes. 

2. VELOCITY OF THE LOCAL STANDARD OF REST 

Despite decades of effort, the local MW rotation rate 
remains poorly know n, with measurements va rying by 
25%. Hipparcos fi's l|Feast fc Whitelocklli997ft suggest 
that the Local Standard of Rest (LSR) velocity is ©lsr = 
(217.5 ± 7.0) (i? /8) km s" 1 — i.e. near the IAU adopted 
value of 220 km s _1 . But a more recent measurement of 
jU for Sgr A* (Reid & Brunthaler 2004) yields a higher 
(235.6 ± 1.2)(i? /8) km s" 1 , whereas direct HST mea- 
surements of the /it's of bulge stars against background 
galaxies in the same field yield (202.4 ± 20.8)(i? /8) km 
s- 1 (Kalirai et al. 2004) and (220.8±13.6)( J R /8) km s" 1 
(Bcdin et al. 2003). Of course, these measures (as well 
as any of those depending on the Oort constants) rely on 
an accurate measure of Rq (§1). The solar peculiar mo- 
tion must also be known, but is a smaller correction (e.g., 
5.3 ± 0.6 km s -1 ; Dehnen & Binney 1998). On the other 
hand, considerations of non-axisymmetry of the disk 
yield corrections to the measurement s that suggest Olsr 
may be as low as 184 ± 8 km s~ x (Oil ing fe Merrifieldl 
1998) or lower (Kuijken & Tremaine 1994). Independent 
methods to ascertain Olsr are of great value because it 
is fundamental to establishing the MW mass scale. 

Eventually, as part of a Key Project of SIM, Olsr 
will be measured directly by the absolute /i of stars near 
the GC. Here we describe an independent method for 
ascertaining Olsr using halo tidal streams that over- 
comes several difficulties with working in the highly dust- 
obscured, crowded GC, and one also insensitive to Rq 
(for all reasonable values of the latter). The ideal tidal 
stream for this method is one with an orbital pole ly- 
ing near the Yqc axis. The Sgr tidal stellar stream not 
only fulfills this requirement, but its stars, particularly its 
trailing arm M giants, are ideally placed for uncrowded 
field astrometry at high MW latitudes, and at relatively 
bright magnitudes for, and requiring only the most mod- 
est precisions from, SIM. Indeed, as we show, this method 
is even within the grasp of future high quality, ground- 
based astrometric studies. 

It is remarkable that the Sun presently lies within a 
kiloparsec of the Sgr debris plane (MSWO). The pole of 
the plane, (l p ,b p ) = (272,-12)°, means that the line of 
nodes of its intersection with the MW plane is almost 
coincident with the Xqc axis. Thus (Fig. 2) the motions 
of Sgr stars within this plane are almost entirely con- 
tained in their Galactic U and W velocity components, 
whereas the V motions of stars in the Sgr tidal tails al- 
most entirely reflect solar motion. To the degree that its 

V distribution is not completely flat in Figure 2 is due to 
the slight amount of streaming motion projected onto the 

V motions from the 2° Sgr orbital plane tilt from Xqc, 
compounded by (1) Keplerian variations in the space ve- 
locity of stars as a function of orbital phase, as well as (2) 
precessional effects that lead to A Q -variable departures 
of Sgr debris from the nominal best fit plane to all of 
the debris. The latter is negligible for trailing debris but 
is much larger for the leading debris, which is on aver- 
age closer to the GC and dynamically older compared to 
the trailing debris when viewed near the Galactic poles 
(J05). In addition, because the leading debris gets arbi- 
trarily close to the Sun (L05), projection effects make it 
more complicated to use for the present purposes. Addi- 
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tional problems with the leading arm debris, which sug- 
gest that more complicated effects have perturbed it are 
discussed in L05 and J05. In contrast, the Sgr trailing 
tail is beautifully positioned fairly equidistantly from us 
for a substantial fraction of its stretch across the South- 
ern MW hemisphere (MSWO). This band of stars arc- 
ing almost directly "beneath" us within the Xqc-Zgc 
plane provides a remarkable zero-point reference against 
which to make direct measurement of the solar motion al- 
most completely independent of the Sun's distance from 
the GC. The extensive mapping (MSWO) of the Sgr 
tails with Two Micron All-Sky Survey (2MASS) M gi- 
ants provides an ideal source list for individual stellar 
targets from this > 360°-wrapped, MW polar ring. 

L05 used M giant spatial (MSWO) and RV data (Ma- 
jewski et al. 2004) to constrain models of Sgr disrup- 
tion, best fitting when a 3.5 xlO 8 M Q Sgr of 328 km s _1 
space velocity orbits with period 0.85 Gyr and apo:peri- 
Galactica of 57:13 kpc. These models fit what appears 
to be ~2.5 orbits (2.0 Gyr) of Sgr mass loss in M giants. 
The adopted MW potential is smooth, static and given 
by the sum of a disk, spheroid and halo described by the 
axisymmetric function <&haio = w haio m (^ 2 + [z / q] 2 + d 2 ) 
where q is the halo flattening, R and z are cylindrical 
coordinates and d is a softening parameter. Additional 
model details are given in L05. 

Model fits to the Sgr spatial and velocity data allow 
predictions of the 6-D phase space configuration of Sgr 
debris. Figure 2 shows predicted U, V, W velocity com- 
ponents of debris as a function of longitude, A©, in the 
Sgr orbital plane (see MSWO). The debris is shown as- 
suming q = 0.9, Ro = 7 kpc, and a characterization of 
the total potential whereby the LSR speed is Olsr = 220 
km s _1 . L05 explores how variations in q affect primarily 
the U and W (through projection along RV). Figure 3 
(green, yellow, and magenta points) shows how variations 
in the scale of the potential, expressed through variations 
in adopted ©lsr, affect V . Clearly, Olsr ranging from 
180 to 260 km s _1 translates to obvious variations in 
observed V for trailing arm stars. This effect is easily 
separable from any residual uncertainty in the shape of 
the potential or i?o : Figure 3 (red and blue points re- 
spectively) illustrates negligible V changes produced by 
holding Olsr fixed at 220 km s _1 but varying q from 0.9 
to 1.25 (i.e. oblate to prolate) and Rq from 7 to 9 kpc. 

Figure 3 is the basis for the proposed use of Sgr to mea- 
sure Olsr- Ideally, to execute the experiment requires 
obtaining V from the observed /j, and RVs of Sgr arm 
stars. However, because of the particular configuration 
of Sgr trailing arm debris, almost all of V is reflected in 
the jU of these stars, and, more specifically, the reflex solar 
motion is contained almost entirely in the /i/ cos(6) com- 
ponent of \x for Sgr trailing arm stars away from the MW 
pole. Working in the observational, [i regime means that 
vagaries in the derivation of individual star distances can 
be removed from the problem, as long as the system is 
modeled with a proper mean distance for the Sgr stream 
as a function of A© . Figure 4 shows three general regimes 
of the trailing arm hi cos(6) trend: (1) A© > 100° where 
Hi cos(6) is positive and roughly constant, (2) the region 
from 100° > A© > 60° where m cos(6) flips sign as the 
debris passes through the South Galactic Pole to shift the 
Galactic longitudes of the trailing arm by <~ 180°, and 
(3) A© < 60°, where /i/cos(fe)is negative and becomes 



smaller with decreasing A© (because the Sgr stream be- 
comes increasingly farther). The sign flip in /i; cos(6) is 
a useful happenstance in the case where one has /z data 
not tied to an absolute reference frame but which is at 
least robust to systematic zonal errors: In this case the 
peak to peak amplitude of /i; cos(6) for the trailing arm 
stars yields (two times) the reflex motion of the Sun 5 . 

The intrinsic RV dispersion of the Sgr trailing arm has 
been measured to be ~10 km s _1 (Majewski et al. 2004); 
assuming symmetry in the two transverse dimensions of 
the stream gives an intrinsic /i dispersion of the Sgr trail- 
ing arm of ~0.1 mas yr _1 (see Fig. 4a). Thus, until SIM- 
quality proper motions exist, the measurement of the re- 
flex solar motion by this method will be dominated by 
the error in /i. To quantify the accuracy of the proposed 
method, we introduce artificial random errors into the 
proper motions of the five models shown in Figure 3 and 
calculate the accuracy with which we expect to recover 
the solar reflex motion. 

Simply applying the formalism described above, we re- 
cover 

Olsr values of 6 212, 255, and 279 km s" 1 m mod- 
els for which input Olsr = 180, 220, and 260 km s _1 
respectively. This indicates that the method systemati- 
cally overpredicts Olsr by about 30 km s _1 ; this is be- 
cause (see Fig. 3) the trend of V with A© is not perfectly 
flat but changes by ~ 30 km s _1 between the peaks at 
A© = 60 — 65° and 115-120°. Correcting for this system- 
atic bias, we perform 1000 tests where we randomly draw 
particles from the model debris streams in these ranges 
with artificially added random scatter in the it, and find 
that recovering the solar velocity to within 10 km s _1 re- 
quires a sample of approximately 200 stars with \x mea- 
sured to about 1 mas yr _1 precision with no zonal sys- 
tematics. Using the models with Olsr = 180/220/260 
km s _1 , these tests recover mean values of 182/225/249 
km s _1 respectively with a dispersion of results between 
the tests of 10 km s . As expected from Figure 3, vary- 
ing q and Rq has negligible effect: Tests on models in 
both of these MW potentials (where Olsr = 220 km 
s _1 ) recover mean values of 225 and 228 km s _1 . 

Present astrometric catalogs are just short of being 
able to do this experiment: Hipparcos is not deep enough, 
the Southern Proper Motion Survey (Girard et al. 2004) 
has not yet covered enough appropriate sky area, and 
UCAC2 (Zacharias et al. 2004) has several times 
larger random errors than useful as well as comparably- 
sized zonal systematic errors at relevant magnitudes (N. 
Zacharias, private communication). However, to demon- 
strate how only modest advances in all-sky jj, precisions 
are needed to make a definitive measurement, Figure 4 
includes a direct comparison of the m cos(6) trend for Sgr 
M giants using UCAC2 /z's for 2MASS M giants. Impres- 
sively, the overall expected fii cos(6) trends can be seen, 
but the large scatter and systematic shifts in the trail- 
ing arm motions belie the limits of UCAC2 accuracies at 
V ~ 15. Even a factor of two improvement in UCAC2 
random errors and elimination of zonal errors might lead 
to a useful measurement of Olsr- It is not unreasonable 
to expect advances in all-sky ji catalogues at this level 

5 We find that these peaks lie at A = 60 - 65° and 115 - 120°; 
note that technically V oc (n/d)eo—65 + (mAOhs— 120 ■ 

6 Correcting for the assumed 12 km s -1 speed of the Sun with 
respect to the LSR. 
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soon (e.g., from the Origins Billion Star Survey or Gaia), 
but in any case SIM PlanetQuest will easily obtain the 
necessary /x (and parallaxes) of selected Sgr trailing arm 
giants. 

We appreciate funding by NASA/JPL through the 



Taking Measure of the MW Key Project for SIM Plan- 
etQuest, NSF grant AST-0307851, the Packard Foun- 
dation, and the F.H. Levinson Fund of the Peninsular 
Community Foundation. SRM appreciates the hospital- 
ity of the Carnegie Observatories during the writing of 
this paper. 



REFERENCES 



Baldwin, J. E., Lynden-Bell, D., & Sancisi, R. 1980, MNRAS, 193, 
313 

Bedin, L. R., Piotto, G., King, I. R., & Anderson, J. 2003, AJ, 126, 
247 

Dehnen, W., & Binney, J. J. 1998, MNRAS, 298, 387 
Eisenhauer, F., Schodel, R., Genzel, R., Ott, T., Tecza, M., Abuter, 

R., Eckart, A., & Alexander, T. 2003, ApJ, 597, L121 
Feast, M. W. & Whitelock, P. 1997, MNRAS, 291, 683 
Gardiner, L. T., & Noguchi, M. 1996, MNRAS, 278, 191 
Girard, T. M., Dinescu, D. I., van Altena, W. F., Platais, I., Monet, 

D. G., & Lopez, C. E. 2004, AJ, 127, 3060 
Johnston, K. V., Law, D. R., & Majewski, S. R. 2005, ApJ, 619, 

800 (J05) 

Kalirai, J. S., et al. 2004, ApJ, 601, 277 
Kuijkcn, K., & Tremaine, S. 1994, ApJ, 421, 178 
Law, D. R., Johnston, K. V., & Majewski, S. R. 2005, ApJ, 619, 
807 (L05) 

Majewski, S. R. 2004, Pub.Astr.Soc. Australia, 21, 197 



Majewski, S. R., et al. 2004, AJ, 128, 245 

Majewski, S. R., Skrutskie, M. F., Weinberg, M. D., & Ostheimer, 

J. C. 2003, ApJ, 599, 1082 (MSWO) 
van der Marel, R. P., Alves, D. R., Hardy, E., & SuntzefT, N. B. 

2002, AJ, 124, 2639 
Newberg, H. J., et al. 2002, ApJ, 569, 245 
Oiling, R. P. & Merrifield, M. R. 1998, MNRAS, 297, 943 
Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 
Richter, O.-G., & Sancisi, R. 1994, A&A, 290, L9 
Rix, H., & Zaritsky, D. 1995, ApJ, 447, 82 
Shapley, H. 1918, ApJ, 48, 154 
Vivas, A. K., et al. 2001, ApJ, 554, L33 

Walker, I. R., Mihos, C, & Hernquist, L. 1996, ApJ, 460, 121 
Zacharias, N., Rafferty, T. J., & Zacharias, M. I. 2000, ASP Conf. 
Ser., 216, 427 



Fundamental Galactic Parameters 



5 





i i 


i i I i i i i 


i i i i i i i i i i 


I i i i i I i 


1 1 1 


i i i 


1.5 










• io 5 


stars — 












■ io 4 


stars - 


1 










▲ io 3 


stars - 








o 


A 






0.5 




Q 


□ 












A 


n 


1 













3 












A 


§ 


O 




0.5 




▲ 


• 














A 








-1 










□ 




1.5 










I 






~ i i 


1 1 i i i i i 


1 i i i i 1 i i i i 


1 i i i i 1 i 


, , i 


i i i 



7 7.5 8 8.5 9 9.5 10 10.5 
log(Progenitor Initial Mass) 



Fig. 1. — Results of plane-fitting to simulations of disrupting satellites of different masses orbiting similarly to the Magellanic Clouds. 
Filled, half-filled and open symbols arc for simulations with no, 10% and 20% distance errors imposed on the member stars. 
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Fig. 2. — Predicted U, V, W velocities (right-handed system) as a function of longitude in the Sgr orbital plane (Aq = 0° at present Sgr 
position). Debris lost on last half (yellow) and previous full (magenta) orbits are shown (see Law et al. 2005). Trailing debris stretches 
from A = 0° to ~ 180° ; leading arm debris goes from A = 0° to ~ 270° . The Galactic model has a q = 0.9 halo with 0lsr = 220 km 
s -1 and Rq = 7 kpc. 
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Fig. 3. — Variation of the Sgr trailing arm V velocities for a range of SlsRi an d 9- Green, yellow, and magenta points represent 
debris from satellites disrupting in potentials where q = 0.9 and Rq = 7 kpc, but 0lsr i s 180, 220, and 260 km s —1 , respectively. Red and 
blue points represent satellite debris in potentials where 0lsr = 220 km s — 1 but q = 1.25 (i?o fixed at 7 kpc) and Rq = 9 kpc (q fixed at 
0.9), respectively. Note how changes in the shape of the potential have little to no effect on V, while changes in the scale of the potential 
produce large, approximately linear shifts. 
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Fig. 4. — (a) Same as Fig. 2 (i.e., same MW model), but for predicted fn cos(fe) and for 2.5 orbits of mass loss, (b) Observed UCAC2 
fil cos(b) for 2MASS M giants in Galactic regions dominated by Sgr stream stars (all M giants within 7 kpc of the nominal Sgr plane and 
having heliocentric distances of 15-30 kpc). Note the differing vertical scales of the two panels. 



